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Abstract 

We find that general relativity can be naturally free of cosmological sin- 
gularities. Several nonsingular models are currently available that either 
assume ad hoc matter contents, or are nonsingular only over a sector of so- 
lution space of zero measure, or depart drastically from general relativity at 
high energies. After much uncertainty over whether cosmological inflation 
could help solve the initial-singularity problem, the prevailing belief today 
is that general relativistic cosmology, with inflation or without, is endemi- 
cally singular. This belief was reinforced by recent singularity theorems that 
take account specifically of inflation. Here, a viable inflationary cosmology 
is worked out that is naturally free of singularities despite the fact that 1) it 
uses only classical general relativity, 2) it assumes only the most generic in- 
flationary matter contents, and 3) it is a theory of the chaotic-inflation type. 
That type of inflation is the most widely accepted today, as it demands the 
least fine-tuning of initial conditions. It is also shown how, by dropping the 
usual simplification of minimal coupling between matter and geometry, the 
null energy condition can be violated and the relevant singularity theorems 
circumvented. 
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The powerful singularity theorems of the sixties [1,2] seemed to leave 
little room for avoiding a catastrophic singularity in the cosmo logical past, 
implying a necessary primeval breakdown of the laws of physics. Naturally, 
questions were then raised about the self-consistency of general relativity as 
a complete theory of gravity, a philosophical question the answer to which is 
still far from clear [3]. 

There have been several successful attempts at building alternative the- 
ories of gravity that are nonsingular in the limit of very high energies (see, 
e.g., [4-6]). It has also been shown that the equations of motion of several 
scalar-tensor theories of gravity may admit nonsingular asymptotic solutions 
if the scalar functions in the Lagrangian are reverse-engeneered for that spe- 
cific purpose [7-9]. Not surprisingly, most of the corresponding full solutions 
are either highly unstable (and hence could not be simulated numerically 
in the literature) or occupy a sector of measure zero in solution space or in 
parameter space. 

Perhaps the most serious prospect of a nonsingular cosmology based on 
ordinary gravity arose from the realization that cosmological inflation [10] is 
probably future-eternal and hence perhaps symmetrically past-eternal as well 
[11-16]. However, this prospect was dimmed considerably by new inflation- 
specific singularity theorems [17-21] that condemned a large class of infla- 
tionary models to be necessarily singular. 

In hindsight, the failure of inflation to guarantee the avoidance of a cos- 
mological singularity is not surprising. The inflationary expansion itself, no 
matter how large, should not affect the global topology of spacetime and 
hence should have no bearing on the existence or non-existence of global sin- 
gularities. However, inflation could potentially play a role in the singularity 
game by the mere fact that it assumes the Universe to be dominated by 
scalar matter at very high energies. Such a universe could easily violate the 
"strong energy condition" which is assumed by several singularity theorems 
[2,22-24]. 

As we shall see below, it is the pre-inflationary dynamics and energetics 
of the scalar matter that are likely to help solve the initial-singularity prob- 
lem. Once it enters the inflationary phase proper ("slow-roll" conditions) 
the scalar field acquires dynamics and energetics that do seem to make the 
model fall under the power of most singularity theorems. More specifically, 
it starts then to obey the "null energy condition." 

We now discuss this latter point more quantitatively. Then, we show in a 
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specific model how this allows an inflationary cosmology to be nonsingular, 
even though the characteristics of its "inflationary stage" proper may suggest 
otherwise. Finally, we follow the evolution of this ordinary-gravity, chaotic in- 
flationary scenario starting from the late inflationary stage, when pre-galactic 
density fluctuations are generated, back through the very brief but crucial 
pre-inflationary phase, and finally through the cosmological bounce itself. 
We confirm numerically that the model is a nonsingular cosmology that is 
stable in solution space. We conclude with some additional comments on the 
relevant recent literature. 

To see how the quasi-de Sitter phase of inflation does not necessarily ef- 
fect the initial-singularity problem, consider first the Raychaudhuri equation 
for null geodesies in the absence of torsion and shear [2]: 

JO 1 

2 + -e 2 = -R, u N»N» . (1) 

9 is the divergence of the considered congruence of null geodesies; v is an 
affine parameter; R^ v is the Ricci curvature tensor; and the iV M are compo- 
nents of a null 4- vector field that is tangent to a Robertson- Walker spacetime 
with metric 

ds 2 = -dt 2 + a 2 {t)dS 2 . (2) 

dS 2 is the metric of a homogeneous and isotropic 3-space. If one chooses the 
null 4- vector field to be locally (1, lt/a(t)), where rf is a unit 3- vector, one 
obtains 

d9 1 „ 9 (a a 2 k\ ( ■ k \ 

where H = a/a is the expansion rate of the spacetime and k = +1 , or — 1 
according to whether the spatial 3-geometry in Eq.(2) is closed, flat or open, 
respectively. Dots denote derivation with respect to t. 

If one assumes general relativity, considers a model that is driven by 
a scalar field <fi with potential V((f)), and makes the usual simplification of 
minimal coupling between <fi and the scalar curvature (£ = in Eq.(lO) 
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below), then (h — k/a 2 ^j is automatically negative definite: 

H-- = -AixGtf , (4) 
or 

where G is the gravitational constant. This seemingly unavoidable fact causes 
enough past convergence through Eq. (3) that past null geodesies are necessar- 
ily incomplete, i.e., the cosmology has a global past singularity [17-21]. Note 
that a de-Sitter spacetime produces zero on the right-hand-side of Eq.(4) 
only if it is exactly de-Sitter, in which case it cannot be part of a consistent 
cosmology. 

The positivity of — (h — k/a 2 ^j , which is essentially the sum of the density 
and the pressure of the model, also implies that the null energy condition 

R^N»N V > (5) 

is satisfied, bringing such a cosmology under the wrath of several singularity 
theorems [17-21]. We comment further on this and other energetic aspects 
at the end of this letter and in [25]. 

The key Eq.(4) above derives from a Lagrangian of the form 

C = — — + -d- u ^ - V(6) , (6) 

where R is the scalar curvature. However, it is necessary for the renormaliz- 
ability of the theory that the Lagrangian contain an explicit coupling term 
between the scalar field and the scalar curvature [26-28]. The correct form 
of the Lagrangian is thus 

Note that in our conventions the conformal value of the nonminimal-coupling 
constant £ is —1/6. Field theoretic considerations do not a priori constrain 
the value of £ (For a more complete discussion, see [29].) However, models 
with £ < are generally too pathological to yield viable inflationary cosmolo- 
gies [29,30-34]. In contrast, it was shown some time ago that models with 
£ > 0, and especially those with £ >> 1, lead to well behaved inflationary 
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cosmologies with several attractive features [30-37], such as the easing of the 
observational constraints on the constants that characterize V((f>). 

It can be shown (Eqs. (26,28) of [35]) that during the slow-roll phase of 
chaotic inflation driven by the theory (7), one has 



Hence, for the cosmologically interesting range £ > 0, one has H < 0, just 
as in Eq.(4) where the effect of nonminimal coupling was neglected. Thus, 
from the above discussion of the Raychaudhuri equation, one could presume 
that the nonminimal-coupling adjustment would have no obvious effect on 
the singularity aspect of the model. However, as we shall find soon, the 
2?re-inflationary stage of the model is drastically effected by the nonminimal 
coupling, allowing e.g. H to assume positive values at the earliest stages, 
and eventually leading to a nonsingular cosmological past. 

Let us then analyze the full equations of motion [32-35] without assuming 
slow-roll conditions. The "energy" and "momentum" equations are, respec- 
tively 




8nG 



Itf + V{<j>) - 6£#00 , (9) 



8ttG 



£#00-£0 2 -£00 



(10) 



The matter field evolves according to 



4> + 3#0 - 6£ \H + 2H 2 + 4 U + 



dV{<j>) 
d(j) 



= . 



(11) 



At this point, one cannot limit oneself to looking for de Sitter-like asymp- 
totic solutions and, if there are any, conclude that the model must be nonsin- 
gular. In fact, many of the candidate nonsingular models in the literature do 
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have two exponentially expanding asymptotic solutions, but fail to cross from 
one asymptotic region to the other, or cross over unstable or zero-measure 
trajectories in phase space [8,9,38,25]. In fact, the present model illustrates 
well the caution that must be exercised in inferring the existence or non- 
existence of singularities from either the properties of the inflationary stage 
alone or those of the asymptotic solutions in isolation. In the present case, 
for example, H < during asymptotic slow-roll phases and yet the model 
can manage to avoid crushing into a singularity by very briefly shifting away 
from the slow-roll conditions, allowing H to turn positive. Again, there 
clearly is not a causal relationship between the existence of de-Sitter like 
asymptotic solutions and the avoidance of cosmological singularities. This 
aspect is not always recognized, although it had become clear during, e.g., 
the investigation of "pre-big-bang" cosmology [38-40,5]. 

Let us choose the potential to be of the simple and generic form V{4>) = 
X(f> 4 . We start at the late stage of inflation when the slow-roll conditions 
hold. Then, the model evolves according to 

H 2 ~ -B 1 log a + d , 4> 2 ~ -B 2 log a + C 2 , a ~ exp (c 3 t - ^-Bi* 2 ) , 

(12) 

where C±, C 2 and C3 are constants that are determined by the initial con- 
ditions, while B\ and B 2 are positive constants that depend on the coupling 
parameters A and £: 

B 2 = ^- , B 1 = —B 2 . (13) 

In addition, <fi ~ —B 2 H/2<p and fluctuations are produced that will eventu- 
ally result (when the relevant scales cross back into the Hubble radius after 
the end of inflation) in pre-galactic density perturbations with magnitude 
[35,36] 

— oc -. . (14) 

p 0(1 + 6O 1/2 
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As announced above (see Eq.(8)), H is negative-definite during this rela- 
tively late stage of inflation: H —Si/2. As one evolves the model towards 
the past, the slow-roll solution (12) breaks down and, a priori, H can cover a 
wide range of values of either sign (see Eqs.(9-10)), depending on the choice 
of parameters and initial conditions. In a separate paper, we explore the 
many possible detailed behaviors of the model that can result [25]. Our fo- 
cus here is the possibility of obtaining a complete and naturally nonsingular 
inflationary model from fairly generic initial conditions. (By this, we mean 
of course generic inflationary initial conditions. The question of whether the 
latter are themselves probable among all possible initial conditions is still 
open [43]). 

By inspecting Eqs.(9-10) analytically for the behavior of the matter func- 
tions (0, 0, 0) and the metric functions (a, H, H), we find that three types 
of scenarii are possible. 

1) H remains negative-definite throughout and H diverges at a finite 
past value of t where a(t) plunges to zero. I.e., one hits a singularity. This 
behavior is generic in the case of an open spatial geometry (k = —1). 

2) H first becomes positive (again, as we run the model toward the past) 
but then decreases toward zero asymptotically. In this case, which obtains 
typically for spatially flat geometries (k — 0), a(t) decreases asymptotically 
toward a finite minimal size. Such solutions are technically nonsingular, but 
are usually unstable for most combinations of parameters. Moreover, in the 
real Universe, the condition of rigorous spatial flatness has arguably zero 
measure, so that unlike in lower-energy cosmological problems, the cases of 
practical interest here are probably only k = — 1 and k — +1. 

3) H eventually becomes positive, and stays so sufficiently far in the past 
to allow H to reach zero and turn negative. The scale factor a(t) reaches a 
minimum and re-expands again toward the past. This cosmological bounce 
obtains for the spatially closed case (k = +1). 

In Figs. (1-3), we confirm numerically the above analysis, focusing on the 
latter nonsingular case. (More output from the first two cases can be found 
in [25].) We simulate that scenario starting from the late inflationary stage, 
back through the change of signs in H, and finally, safely through the cos- 
mological bounce itself. Fig.(l) shows the evolution of the scalar field in 
the neighborhood of the bounce. Fig. (2) displays the overall behavior of 
the expansion rate H, including the inflationary phases (straight portions of 
the curve, see Eq.(12)) and the non- inflationary, or inter-inflationary phases 
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where H turns positive and the null energy condition is violated. Fig. (3) 
shows clearly the distinction between inflationary and inter-inflationary evo- 
lutions of the scale factor, and also confirms that the model produces enough 
expansion to be observationally viable. To facilitate the comparison with the 
published literature, the figures shown here are obtained for £ = 1000 and 
A = 0.001, just as in [32,33,35]. Qualitatively, the results hold for a broad 
range of parameter values. One can also obtain some potentially interesting 
variations by adding a mass term to the Lagrangian or by specializing to 
very small or very large parameter values [25]. Here, our aim is to show that 
a generic case can be nonsingular. 

We have also checked the special case of minimal coupling (£ = 0). One 
major difference in this case is that the null energy condition can never be 
violated (see Eqs.(4,5)), so that most of the singularity theorems do apply. 
The sub-cases k — and k = — 1 are always singular, as the theorems 
predict [17-21]. For k = +1, the situation is less trivial, because technically 
the equations of motion do not exclude a bounce. (Note that the inflation- 
specific singularity theorems are also more technically involved in this case.) 
For instance, one can set H = <fi = with £ = and k — 1 in Eqs.(9-ll), 
and successfully evolve the system out of that state of minimal size in both 
directions of time. Unfortunately, because of Eq.(4), H turns negative almost 
immediately, and the model recollapses far too soon to constitute a viable 
inflationary scenario. One could still obtain a longer period of expansion by 
forcing to remain virtually zero for a long enough time. For example, this 
could be arranged by using a potential that is extremely flat near = or 
that has a local minimum around that value, and placing the initial <fi and 
very close to zero. But one would then be resorting to an inflation of the 
"new" or the "old" type, respectively. Compared to chaotic inflation, those 
two types have many additional observational and conceptual difficulties and 
are considered today to be among the least attractive [10,41-43]. 

Still on the topic of minimal coupling, a very recent work was just brought 
to our attention which explores the energy constraints that an inflationary 
model must satisfy in order to allow for a cosmological bounce [44]. Unfor- 
tunately, that study misses the possibility (probably the necessity [26-28]) 
of nonminimal coupling (e.g., stating inaccurately that all inflationary mod- 
els are minimally coupled.) The study also concludes that only "old" and 
"new" inflation ("and not chaotic or eternal inflation") are compatible with 
the bounce scenario from an energetics point-of-view, which is also inaccurate 
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as we saw above. Finally, reference [44] does not actually work out a specific 
inflationary model. Doing so would be necessary because, as we indicated 
above, one can easily construct solutions that look like they could bounce, ei- 
ther because of their asymptotic behavior or on energetics grounds, but that 
are not viable cosmologies [7-9,38,25]. Still, [44] is an interesting extension 
of previous investigations of the interplay between the energy conditions, the 
singularity theorems, and inflation. 

Our results suggest the possibility of a Universe that inflates, then re- 
heats (when drops close enough to zero) and expands as a Friedmann 
universe, then recontracts and reheats again, becomes eventually scalar-field 
dominated, then bounces when it reaches a relatively large minimum size, 
then inflates and reheats again, and so on for ever in both directions of 
time (see Figs. (2,3)). Of course, one would have to contend then with the 
rapid growth of anisotropies and inhomogeneities during contracting phases, 
and with whether black hole formation can be kept in check, e.g., by the 
possible dominance of some type of dark matter, in particular some weakly- 
interacting, nonbaryonic dark matter. Alternatively, it is perhaps possible 
to cast these findings in the context of "eternal" inflation [11-16], and ask 
whether the latter can thus be made eternal in the past as well as in the 
future [21]. 

We hope to have shown that it might not be necessary to depart much 
from general relativity to avoid cosmological singularities, were they to be 
eventually deemed physically or philosophically unacceptable. If, as it is 
suggested here, the minimal size to which the Universe recollapses is not too 
small, it might not be necessary to call upon quantum gravity in order to 
build fully consistent cosmologies. 

During the course of this research, I have benefited from many instructive 
exchanges on related technical and philosophical aspects with R. Branden- 
berger, V. Husain, E. Luft, S. Savitt, K. Schleich, W. Unruh, P. Varniere, A. 
Vilenkin, S. Weinstein, and D. Witt. 
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Figure 1 
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Figure 2 
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Figure 3 
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